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A  numerical  investigation  of  high-resolution  schemes  for  solving  the  compressible  Navier-Stokes 
equations  in  the  context  of  Implicit  Large  Eddy  Simulation  (ILES)  is  presented.  Computations  have  been 
performed  using  a  hybrid  total  variation-diminishing  (TVD)  scheme.  The  hybrid  TVD  scheme  combines 
third-order  Godunov-type  fluxes  with  first-order  (dissipative)  schemes  using  the  flux  limiters  approach. 
Results  are  presented  from  numerical  experiments  of  turbulent  flow  past  an  open  cavity. 


1  INTRODUCTION 

It  has  been  observed  for  more  than  a  decade  that  high-resolution  methods  have  characteristics  that  mimic 
the  effects  of  finite  viscosity  and  thermal  conductivity  and  lead  to  a  ‘dissipation’  which  is  similar  to  that 
achieved  by  adequately  resolved  numerics  involving  finite  Reynolds  and  Prandtl  numbers  [1],  [2],  [3],  [4], 
[5],  [6],  The  idea  to  use  these  methods  as  an  implicit  way  to  numerically  model  complex  turbulent  flows, 
e.g.,  flows  dominated  by  vorticity  leading  to  turbulence,  flows  featuring  shock  waves  and  turbulence,  and 
the  mixing  of  materials,  is  an  evolving  area  of  research  referred  to  as  Implicit  Large  Eddy  Simulation 
(ILES)  or  Monotonically  Integrated  LES  (MILES).  Such  flows  are  extremely  difficult  to  practically  obtain 
stably  and  accurately  in  under-resolved  conditions  (with  respect  to  grid  resolution)  using  classical  linear, 
both  second  and  higher-order  accurate  schemes. 

The  success  of  high-resolution  methods  to  compute  turbulent  flows  seems  to  depend  on  a  delicate  balance 
of  truncation  errors  due  to  wave-speed-dependent  terms  (chiefly  responsible  for  numerical  dissipation)  in 
the  case  of  Godunov-type  methods  and  hyperbolic  part  of  the  flux.  It  is  the  essence  of  this  balance  that 
needs  to  be  understood  in  order  to  regularise  the  flow  and  allow  shock  propagation  to  be  captured 
physically  realistically  even  if  it  is  not  fully  resolved  on  the  computational  mesh.  Nonlinear  limiters  in 
high-resolution  methods  guard  the  methods  against  catastrophic  failures  due  to  nonlinear  wave-steepening 
or  unresolved  features  by  triggering  entropy  producing  mechanisms  that  safeguard  the  calculation  when 
the  need  arises.  The  two  key  questions  are:  (i)  what  criteria  should  be  used  to  design  the  nonlinear  limiter 
that  triggers  the  entropy  production,  and  (ii)  to  what  extent  numerical  dissipation  accounts  for  turbulent 
flow  effects. 

Results  from  the  implementation  of  high-resolution  methods  in  wall-bounded  flows  [3],  [6]  show  that  in 
principle  there  is  nothing  that  prevents  the  use  of  these  methods  in  near  wall  flows  even  without  using  an 
explicit  turbulence  model.  Both  ILES  and  traditional  LES,  based  on  an  explicit  subgrid  scale  (SGS) 
turbulence  modelling,  pose  substantial  challenges  in  high-Reynolds,  near-wall  flows,  especially  in  the 
presence  of  separation  from  gently  curved  surfaces,  where  resolution  and  thus  computing-cost  issues  are 
critical. 

In  the  present  study,  we  investigate  the  ILES  approach  in  the  context  of  under-resolved  simulations,  using 
a  hybrid  TVD  scheme  that  combines  third-order  Godunov-type  fluxes  with  first-order  (dissipative)  fluxes 
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ORGANIZATION 


through  the  implementation  of  limiters.  Results  are  presented  from  a  numerical  experiment  of  turbulent 
flow  past  an  open  cavity. 

2  NUMERICAL  FRAMEWORK 

The  physics  of  (Newtonian)  fluid  flow  is  governed  by  the  Navier-Stokes  equations.  These  equations 
consist  of  the  continuity,  momentum  and  energy  laws  which  can  be  written  in  a  coupled  generalised 
conservation  form 


—  +  V-(/9U)  =  0, 

dt 

+  V  •  (yOUll)  =  -V  •  P, 
dt 

de 

—  +  V  •  (eu)  =  -V  ■  (u-  P)  -  V  ■  q, 
dt 

where  u ,  p  ,  e ,  and  q  stand  for  the  velocity  components,  density,  total  energy  per  unit  volume,  and  heat 
flux,  respectively.  The  tensor  P  for  a  Newtonian  fluid  is  defined  by 

P  =  V{p,T)\  +  |//(V  •  u)I  -  //[(V u)  +  (Vu)T]  , 

where  p (p,  T )  is  the  scalar  pressure,  I  is  the  identity  tensor,  T  is  the  temperature,  and  p  is  the  dynamic 
viscosity  coefficient.  The  above  system  is  completed  by  an  equation  of  state.  For  a  perfect  gas  the  equation 
of  state  is  given  by  p  =  pRT ,  where  R  is  the  gas  constant.  For  the  solution  of  the  Navier-Stokes 
equations,  we  have  employed  an  explicit,  third-order  TVD  Runge-Kutta  scheme  [7],  central  differences 
for  the  viscous  terms  and  a  high-resolution  scheme  for  the  discretisation  of  the  advective  terms.  This 
scheme  is  briefly  described  below  (for  further  details  see  [8],  [9])  for  the  one-dimensional,  inviscid 
counterpart  of  the  equations  in  matrix  form 


au  SE 

dt  dx 


where  U  is  the  array  of  the  unknown  variables  and  E  is  the  flux  associated  with  the  terms  in  x  - 
direction.  The  advective  flux  derivative  5E  /  dx  (similarly,  for  the  other  advective  flux  derivatives)  is 
discretised  at  the  centre  of  the  control  volume  (i)  using  the  values  of  the  inter-cell  fluxes,  i.e., 

dE  _  E.+1/2  —  E._1/2 
dx  Ax 


The  high-resolution  method  employed  here  for  the  calculation  of  the  convective  fluxes  is  a  hybrid  TVD 
scheme.  This  calculates  the  flux  as 


^i+l/2  W  i+ 1/2® 


SW-FVS 

1+1/2 


+  (1-^/+l/2)E/ 


CB 

+1/2 


where  E^v[iS  and  E^/2  are  the  inter-cell  fluxes  according  to  a  Steger- Warming  flux  vector  splitting 
(SW-FVS)  scheme  and  a  Godunov-type,  characteristics-based  (CB)  scheme,  respectively.  The  term  i//I  ll/2 
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is  a  limiter  function  [8]  defined  by  the  square  of  the  (local)  Mach  number  differences  across  cell  faces. 
Second-  and  third-order  of  accuracy  is  achieved  in  discontinuities  and  smooth  flow  regions,  respectively. 

The  Godunov-type,  characteristics-based  scheme  [10]  uses  up  to  third-order  (non-oscillatory)  interpolation 
[8]  to  compute  the  conservative  variables  along  the  characteristics.  The  scheme  is  third-order  accurate  in 
both  space  and  time  when  it  is  used  in  conjunction  with  a  third-order  TVD  Runge-Kutta  scheme  for  the 
time  integration. 

The  SW-FVS  scheme  defines  the  inter-cell  flux  as 

EfXfra=Er+1/2(Ui)  +  E:+1/2(us)  , 

where  the  left,  U  L ,  and  right,  U  R ,  states  of  the  conservative  variables  can  be  obtained  by  first-,  second- 
or  higher-order  accurate  interpolation  schemes.  Here,  the  positive  and  negative  fluxes  are  computed  by  an 
improved  version  of  the  Steger- Warming  flux  vector  splitting  (SW-FVS)  scheme  [8],  [9]  in  conjunction 
with  first-order  interpolation  for  U  L  and  U  R .  The  first  order-interpolation  was  selected  in  order  to 
provide  a  lower-order  (first- order)  flux  in  the  TVD  scheme. 

Limiters  are  the  general  nonlinear  mechanism  that  distinguishes  modem  methods  from  classical  linear 
schemes.  Their  role  is  to  act  as  a  nonlinear  switch  between  more  than  one  underlying  linear  methods  thus 
adapting  the  choice  of  numerical  method  based  upon  the  behaviour  of  the  local  solution.  Limiters  result  in 
nonlinear  methods  even  for  linear  equations  in  order  to  achieve  second-order  accuracy  simultaneously 
with  monotonicity.  Numerical  flux  limiters  can  act  like  dynamic,  self-adjusting  models,  modifying  the 
numerical  viscosity  to  produce  a  nonlinear  eddy  viscosity  [4],  [6]. 


3  RESULTS 

Compressible,  turbulent  flow  past  an  open  cavity  features  a  variety  of  flow  phenomena  including  large  and 
small  vortical  structures,  free  shear  layers,  transitional  flow,  flow  separation  and  flow  re-laminarisation, 
shock  and  rarefaction  waves. 

Under-resolved  simulations  at  Mach  number  0.8  and  Reynolds  number  2500  have  been  carried  out, 
showing  the  development  of  an  oscillating  unsteady  pressure  field,  see  Figures  1-3.  It  is  self-sustained  and 
coherent  but  shifted  in  phase  between  the  cavity's  leading  edge  and  trailing  edge.  The  arising  acoustic 
radiation  is  shown  in  Figure  4.  Acoustic  waves  emanating  from  the  trailing  edge  propagate  upstream  and 
trigger  vortex  shedding.  This  feedback  mechanism  is  responsible  for  the  self-sustained  oscillating  pressure 
field. 

Furthermore,  the  ILES  results  for  the  dominant  Strouhal  number  were  found  in  good  agreement  with  DNS 
of  [11]  for  computations  at  different  Mach  number  conditions,  e.g.,  at  Mach  number  0.8  ILES  and  DNS 
predicted  values  of  0.635  and  0.65,  respectively,  while  at  Mach  number  0.6  the  values  of  the  Strouhal 
number  were  computed  as  0.683  and  0.7,  respectively.  The  comparison  of  dominant  frequencies  shows 
the  applicability  of  ILES  for  open  cavity  flows.  Simulations  for  the  transonic  cases  using  a  SGS 
(Smagorinsky-type)  model  have  shown  no  further  improvement  of  the  results.  This  agrees  with  the 
conclusion  of  [5]  that  when  no  explicit  subgrid  scale  modelling  is  added,  a  high-resolution  method  will  not 
add  any  unnecessary  diffusion. 
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Figure  1:  Pressure  at  the  upstream  face  of  the  cavity. 


Figure  2:  Pressure  at  the  bottom  of  the  cavity. 
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Y=0.005 

Y=0.995 


Figure  3:  Pressure  at  the  downstream  face  of  the  cavity  and  in  the  shear  layer  (X=1.96,  Y=1.00). 


Figure  4:  Instantaneous  density  contours  in  the  computational  domain. 
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SYMPOSIA  DISCUSSION 

REFERENCE  AND/OR  TITLE  OF  THE  PAPER:  34 

DISCUSSOR’S  NAME:  E.  Dowell 
AUTHOR’S  NAME:  M.  Hahn 

QUESTION: 

If  the  grid  is  refined  indefinitely  does  the  LES  solution  approach  the  exact  solution  of  the  Navier-Stokes 
equations? 

AUTHOR’S  REPLY: 

Not  necessarily,  we  still  have  a  certain  amount  of  numerical  dissipation  through  the  truncation  terms.  The 
general  idea  of  LES  is  however  not  to  replace  DNS  or  go  to  very  fine  grids  since  this  is  not  feasible  in  the 
foreseeable  future  due  to  limited  computing  power. 
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